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® ^9?I2 E r- D !f VALUATION DE LA DEGRADATION D'UNE IMAGE VIDEO INTRODUITE PAR UN SYSTEME DE 
CODAGE ET/OU DE STOCKAGE ET/OU DE TRANSMISSION NUMERIQUE. 



(57) L'invention conceme un procede devaluation de la 
degradation d'une image video cod6e par blocs de points 
images ou pixels. II comporte les etapes suivantes: 

a) selectionner une image d'entree et determiner son ac- 
tivity spatiale (SA) dans une fenetre d'analyse selon les 
sous-etapes suivantes: 

i) determiner pour chaque bloc (n f m) les coefficients 
transformers Fn, m (i, y) par une transformed par blocs 

ii) determiner & partir des coefficients transformes Fn, m 
(i, y) I'activite spatiale bsa de chaque bloc 

iii) determiner I'activite spatiale globale SA 1 de la fenetre 
d'analyse 

b) selectionner I'image de sortie correspondant & I'ima- 
ge d'entree et determiner ladite activite spatiale globale SA 2 
de rimage de sortie 

c) comparer les activites spatiales (SAg) et (SA^. 
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PROCEDE D * EVALUATION DE LA DEGRADATION D 1 UNE IMAGE VIDEO 
INTRODUITE PAR UN SYSTEME DE CODAGE ET/OU DE STOCKAGE 
ET/OU DE TRANSMISSION NUMERIQUE 

i^j-^o^hi-^ xuvciiuiuii a uuui uu]tJL uii proceae 
d' evaluation de la degradation d'une image video, 
introduite par un systeme de codage et/ou de stockage 
et/ou de transmission numerique, notamment un systeme 
mettant en oeuvre un codage et/ou un stockage et/ou une 
transmission de signal video numerique avec un bas debit. 

Les degradations sont en general dues aux 
erreurs introduites par les algorithmes de reduction de 
debit et/ou par les liaisons defectueuses de transmission 
ou bien encore aux defauts des codeurs et des decodeurs. 

La numerisation des signaux video ,a ouvert la 
possibility de pouvoir copier, stocker ou transmettre 
1 1 information video en maintenant une qualite constante. 

Cependant, la grande quantite d ' information 
vehiculee par les images video necessite en pratique 
1 1 utilisation de methodes de compression numerique pour 
redduire le debit binaire. 

Une methode de compression tres repandue dans 
le domaine de la video est decrite dans la norme ISO/CEI 
13918 MPEG2. Cet algorithme est dit du type "avec 
pertes", car 1 1 image restituee apres le decodage n'est 
pas identique a 1' original. Afin de maintenir une qualite 
acceptable pour le telespectateur final, les algorithmes 
de reduction de debit tiennent compte des proprietes 
perceptuelles du systeme humain de vision. En depit de 
ceci, les contraintes imposees de debit ou de largeur de 
bande disponible pour la transmission ainsi que le 
contenu du signal video impliquent 1' apparition de 
degradations caracteristiques sur le signal apres 
decodage. Ces degradations introduites par la chaine, par 
exempie MPEG2, influent directement sur la qualite perpue 
de l 1 image finale. 
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L ' evaluation automatique de la qualite des 
signaux audiovisuels a un large champ d ' appl icat ions dans 
ia chaine de television numericue : production, 
distribution, et evaluation des performances des 
5 systemes. 

Par axlleurs, les dispositifs existants ont 
ete elabores pour des tests en laboratoires et ne sont 
pas adaptes pour la teiesurveillance des reseaux de 
distribution. 

10 Des sequences d* images codees MPEG et 

transmises a bas debit pour la diffusion de la television 
numerique ou pour d'autres applications multimedia, vont 
presenter un certain nombre de defauts ou de deformations 
par rapport a la sequence originale. En effet, une liste 
15 non exhaustive de degradations visibles peut etre 
dressee. Les plus perceptibles d'entre elles sont les 
erreurs granulaires, les deformations de contours, les 
pertes d • information, les contours "exot iques " , les 
effet s de bloc, etc.. Cependant, des erreurs de 
20 transmission peu importantes peuvent se traduire par des 
effets plus ou moins localises sur 1 ' image . En cas de 
perturbation importante, elles peuvent se traduire par 
des difficultes d'acces aux informations, par exemple 
coupures de service cu gels d* images pendant un temps 
plus ou moins long selon ia perturbation. L'envergure des 
erreurs depend du niveau de pertinence et de la structure 
des donnees qu * elles affectent : mots de synchronisation, 
vecteurs mouvement, images codees avec ou sans prediction 
ou images de base pour les predictions. En plus des 
coupures ou gels d' images, les degradations constatees se 
materialisent par des blocs ou des macro-blocs errones ou 
mal positionnes. Ceci a pour effet de propager les 
degradations sur toute ia sequence video jusqu'a 1 ' image 
codee sans prediction, done codee independamment des 
35 autres. 



25 



30 



BNSDOCID: <FR 2769453A1 I > 




2769453 



3 

Un procede ci 1 evaluation a ete propose par le 
NTIA (National Telecommunications and Information 
Administration) dans l'article de A. A. WEBSTER et 
Collaborateurs, intitule "An objective videoquality 
5 assesment system based on human perception" et publie en 
juin 1993 dans la revue 5PIE vol.13, p. 15-26. 

Ce procede met en oeuvre une analyse des 
images degradees et des images d'origine, apres les avoir 
filtrees par deux operateurs de SOBEL (matrices 3x3) 
10 vertical et horizontal. L ' image filtree est obtenue par 
convolution en faisant glisser les matrices de SOBEL 
horizontalement et verticalement, et les resultats 
obtenus representent les gradients verticaux et 
horizontaux de 1* image. En d'autres termes, 1 ' image 
15 filtree met en valeur les contours verticaux et 
horizontaux contenus dans 1 1 image initiale (non filtree). 

Une mesure basee sur cette information permet 
de mettre en evidence le changement de contenu entre 
I 1 entree du systeme video et sa sortie. 
20 Le procede propose par le NTIA met en oeuvre 

deux parametres : 

- d'une part, 1 1 inf ormat ion spatiaie SI qui 
represente 1'ecart moyen mesure sur les pixels de 1 1 image 
filtree par i'operateur de SOBEL . II s'agit ici de 
25 determiner l'ecart moyen au niveau des contours de 
1 1 image filtree, en considerant que les contours sont 
importants pour la vision et qu'ils sont affectes par les 
differents traitements des systemes numeriques a 
reduction de debit ; 
30 - d 1 autre part, 1 1 inf ormat ion temporelle TI 

qui represente I'ecart-type de 1 1 image de difference 
entre deux images successives, cet ecart-type etant 
calcule a partir des differences entre les valeurs des 
memes pixels de deux trames successives. Le parametre TI 
35 peut reveler un mouvement saccade du a un defaut du 
codeur . 
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Le procede propose par le NTIA met en oeuvre 
un calcul comparatif cies parametres SI en TI sur le 
signal video numerique, entre une image d 1 entree et une 
image de sortie d'un systeme. 
5 Ce procede preseme un certain nombre 

d' inconvenients . 

Le filtre de SOBEL conserve seulement 
certaines frequences de 1 1 image qui permettent de prendre 
en compte la perte de contour", ce qui fait que la perte 

10 de definition qui pourrait etre prise en compte est celle 
qui se situe dans la gamme des frequences conservees. En 
d'autres termes, la perte de definition ne peut etre 
prise en compte que de maniere partieile. 

D* autre part, le parametre SI prend en compte 

15 des defauts qui tendent a se compenser. En effet, la 
perte d 1 information de 1 1 image tend a faire diminuer le 
parametre SI alors que les faux contours et les effets de 
bloc tendent au contraire a le faire augmenter, ce qui 
fait que le parametre SI n'est significatif que si 1 1 un 

20 ou 1' autre des phenomenes est dominant. 

Enfin, la methode de calcul des parametres SI 
et TI, en utilisant un ecart type calcule de maniere 
globale sur 1 1 image ent iere, reduit de maniere drastique 
1' impact sur ces parametres des degradations localisees. 

- 5 presente invention a pour objet un procede 

qui per met de remedier au moins en partie aux defauts 
precites . 

Le procede selon la presente invention 
exploite une transformee par blocs, par exemple la 

30 transformee en cosinus discrete utiiisee notamment dans 
le standard MPEG, pour met t re en valeur des signatures 
caracteristiques des defauts identifies. Cette approche 
originale permet non seulemen- de mesurer finement les 
erreurs introduites, mais elle permet egalement de tenir 

35 compte du contenu initial du signal video et des 
algorithmes mis en oeuvre dans MPEG. 
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Les transformations par blocs d'une image 
(transformee de Fourier, transformee en cosinus discrete 
TCD, etc..) sont obtenues par 1' operation : 

[Fn,ml = [Tl . [fn,m] . [Tl c ou f(x,y) designe Le 
5 bloc d 1 image a transformer et T(x,y) la matrice de la 
transformation. Une autre transformee par blocs est 
realisee a partir de la transformee en ondelettes de 
1' image, en reorganisant les coefficients d' ondelettes de 
maniere a obtenir des blocs transformes ayant la taille 
10 desiree et en particulier de meme taille que les blocs 
obtenus par les methodes precitees. Un tel mode de 
reorganisation est indique dans l 1 article de 
R. de QUEIROZ et Collaborateurs, intitule "Wavelet 
transforms in a JPEG - like Image Coder", publie en avril 
15 1997 dans la revue IEEE Transactions On Circuits and 
Systems for Video Technology, vol. 7 n°2 p. 419-424. 

L'idee de base de 1' invention est en 
particulier de realiser des calculs sur les blocs selon 
lesquels 1 1 image transmise a ete codee, de maniere a 
20 generer un parametre significatif exempt de l'effet de 
bloc . 

L 1 invention concerne ainsi un procede 
d' evaluation de la degradation d'une image video codee 
par blocs de points images ou pixels, cette degradation 

25 etant ceneree par un systeme de codage et/ou de stockage 
et/ou de transmission produisant une image de sortie a 
partir d'une image d' entree, caracterise en ce qu'il 
comporte les etapes suivantes : 

a) selectionner une image d 1 entree et 

30 determiner une activite spatiale (SA) de 1 1 image d 1 entree 
dans une fenetre d 1 analyse representant au moins une 
partie de 1 ' image qui presente un ensemble de dits blocs 
de points images ou pixels, cette determination mettant 
en oeuvre les sous-etapes suivanres : 

35 i) determiner pour chaque bloc (n,m) du dit 

ensemble de blocs de pixels, les coefficients transformes 
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Fn,m (i,j) par une transformee par blocs selon les dits 
blocs 

ii ) determiner a partir des coefficients 
transformes Fn,m (i,j), l'activite spatiale bsa ce chaque 

5 bloc dudit ensemble de blocs de pixels 

iii ) determiner a partir de l'activite 
spatiale de bsa de chaque bloc l'activite spatiale 
globale SA X de 1 ' ensemble des blocs de pixels constituant 
la fenetre d ! analyse 

10 b) selectionner 1 ' image de sortie 

correspondant a 1 ' image d' entree et determiner ladite 
activite spatiale globale SA 2 de 1 1 image de sortie dans 
ladite fenetre d' analyse, en mettant en oeuvre les sous- 
etapes a i) a a iii) 

15 c) comparer I'activite spatiale globale (SA 2 ) 

de 1 ' image de sortie dans la fenetre d' analyse et 
I'activite spatiale globale (SA X ) de 1 ' image d'entree dans 
la fenetre d f analyse. 

L'activite spatiale bsa d'un bloc (n,ai) peut 

20 etre determinee par la formule suivante : 



25 



35 



f - 



ZZ F ^'A -(^-cao> 



L'activite spatiale bsa peut etre obtenue 
comme indique ci-dessus par combinaison des carres des 
composantes de 1 1 image a partir de la quasi-totaiite des 
composantes de la TCD. Toute autre fonction de 
30 combinaison des composantes de la transformee utilisee 
est applicable pour caracteriser le contenu de 1 1 image . 
Une fonction plus generale est : 



bsa njn = ( ii(k(hj)-F nJ „(i,j))\ -(W).^(0,0))' 

ou k(i,j) est un coefficient constant de ponderation de 
la composante i,j utilisee, et p une constante. 
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II y a plusieurs possibilities pour choisir les 
constantes k(i f j), ez le choix de 1 ' un ou i'autre se fait 
en fonction de 1 1 application recherchee. En effet, la 
fonction k{i,j) est choisie suivanr la caracteristique a 

5 mettre en exergue : tenir compile du systeme visuel ou 
extraire une partie de 1 1 information pertinente. Les cas 
suivants sont proposes : 

1) le procede numero un est a utiliser iorsque 
l'on privilegie une bonne correlation avec la perception 

10 humaine. 

Les valeurs des constantes k(i,j) sont 
initialisers en fonction de 1 1 importance relative des 
coefficients de la transformee pour I'oeil humain, 
notamment sa sensibilite f requentielle, afin de fournir 
15 un parametre d'activite representatif de ce qui est 
pergu. Par exemple, dans le cas de la TCD, on prend 
Jc(i,j) = 1/Q(i,j). Les Q(i,j) sont: les composantes de la 
matrice de quantification exploitee pour la reduction de 
debit, indiquees dans le document suivant qui est extrait 
20 de la norme MPEG-2 : ISO/IEC CD 13818-2 : "Information 
technology - Generic coding of moving pictures and 
associated audio information - Part 2 : video", Dec 1993, 
p 45, § 6.3.7. 

2) le second procede esi utilise iorsque les 
25 coefficients par exemple TCD affectes par la cegradation 
sont identif iables . Les constantes k(i,j) sonz choisies 
de maniere a eliminer certains coefficients de la 
transformee utilisee. En effet, la ponderation k(i,j) est 
utilisee pour retenir les coefficients les plus affectes 
30 ou sensibles a une degradation donnee. II s'agit d'une 
matrice binaire affectant zero aux coefficients a 
eliminer et un aux coefficients pertinents. La selection 
se base soit (a) sur la position du coefficient dans la 
matrice par exemple TCD, soit (b) sur son amplitude 
35 moyenne : 
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a) les coefficients correspondant aux hautes 
frequences spatiales sont souvent ies plus affectes par 
la compression. Un exempie de matrice de ponderation 
suivant le rang du coefficient: par exempie TCD est donne 

5 dans ie tableau ci dessous : 

00001111 
00001111 
0 0 0*0 1111 
10 00001111 

11111111 
11111111 
11111111 
11111111 

15 

b) certains coefficients de faible amplitude 
sont ramenes a zero lors de l'etape de compression. 

Afin de choisir ces coefficients, une moyenne 
ponderee de chacun des coefficients sur la zone d' image 
20 analysee (de taille M.N. blocs) est effectuee : 

.-(• \ ¥ yl ^nj^/} ! , ou Q(i,j) est defini 
m=u=u Q\i 9 j) ! comme ci-dessus, en 1) 

25 Les coefficients dont les moyennes sont parmi 

ies 48 plus faibles vaieurs sont retenus . Pour ces 
derniers, k(i,j)=i et pour les autres k(i,j)=0. 

Dans le cas du parametre "activite spatiale", 
on fixe k(i,j)=l, et p=2 pour la suite du document 
30 relative a la description des figures. 

L'activite spatiale globale SA de i* ensemble 
des blocs peut etre alors determinee par la formule 

suivante : ^ w _ x 

SA = osa „ ^ = V* y 1 bsa _ m 

(HxW) representant le nombre de blocs de 
pixels dans la fenetre d' analyse. 
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Ladite comparaison (etape c) est 

avantageusemenc effectuee a 1'aide du parametre LR defini 
de la maniere suivante : 

LR = gj[fi{SAi,SA2)] avec i,j^{l,2) 
avec /t x-y ou / 2 (x,y)=x/y 

et g l (x)=lOO.|x| ou ft(x)*10a|loB(W)| 
et par exemple : LR = log 10 (SAe/SAs) . 

Le procede peut etre alors caracterise en ce 
qu'il met en oeuvre le calcul du parametre LR sur une 
sequence de M images d' entree et de M images de sorties 
correspondants, et en ce qu'il met en oeuvre le calcul 
d'un parametre d' evaluation ou note de qualite MLR defini 
de la maniere suivante : 

15 MLR = Maximum^ (LR) 

Le procede peut etre avantageusement 
caracterise en ce qu'il comporte la determination d'un 
indicateur d'activite temporelle TA pour un groupe de M 
images, determine de la maniere suivante : 



10 
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f sa<u), pour u variant de 0 a M-l, designant 
ies M coefficients d'une transformee par blocs, par 
exemple une transformer en cosinus discrete appliquee a M 
echantillons temporels de 1'activite spatiale SA. 

Le procede peut etre caracterise en ce qu'il 
comporre, en vue de la determination de l'effet de bloc, 
les etapes suivantes : 

d) determiner l'activite spatiale SA dl de 
1* image d' entree dans une fenetre d' analyse decalee 
presentant des blocs de points image ou pixels qui sont 
decales d'au moins un pixel dans la direction des lignes 
de 1* image et/ou dans une direction perpendiculaire a la 
direction des lignes de 1' image par rapport auxdirs blocs 
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de pixels de ladite fenetre d'analyse, cette 
determination mettant en oeuvre les sous-e:apes suivantes 

i) appliquer ladite transformee par blocs sur 
chaque bloc decale (n,m) dudit ensemble de blocs de 

5 pixels de la fenetre d'analyse decaiee pour determiner 
les coeeficients transformes F d n,m (i,j), \ 

ii) determiner a partir des coefficients 
transformes F d n,m (i,j) des blocs decales, 1'activite 
spatiale bsa d de chaque bloc dudit ensemble de blocs 

10 decales, 

iii) determiner 1'activite spatiale globale 
SA dl de 1» ensemble des blocs decales constituant -la 
fenetre d' analyse decaiee, 

e) determiner 1'activite spatiale globale SA^ 
15 de 1' image de sortie dans ladite fenetre d 1 analyse 

decaiee en mettant en oeuvre pour 1 1 image de sortie les 
sous-etapes d i) a d iii) 

f) comparer d'une part 1'activite spatiale 
globale SA d2 de 1 1 image de sortie dans la fenetre 
d'analyse decaiee et 1'activite spatiale globale SA 2 de 
1' image de sortie dans la fenetre d'analyse, pour evaluer " 
l'effet de bloc dans 1 ' image de sortie et d 1 autre part 
1'activite spatiale SA dl de l 1 image d'enuree dans la' 
fenetre d'analyse decaiee et 1'activite spatiale SA, de ' 
1' image d' entree dans la fenetre d'analyse pour evaluer . 
l'effet de bloc dans 1 ' image d'enrree. 

L'effet de blocs se caracterise de deux I 
manieres differentes : 

intrinsequement, c'est-a-dire par une ; 
fonction representative du contenu de 1 ' image et de 
1 1 image seule. 

SM 2 = gj [/,(SA^SA z )] 
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- de maniere dif f erenciee, c'est-a-dire qu'il. 
est necessaire de comparer deux valeurs d'une fonction 
representative du contenu de 1 1 image : la premiere 
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calculee sur 1 • image de reference (image d'entree) et La 
seconde sur I ' image cegradee (image ce sortie). 

BM = v. \ f ( f (SA .. TJ. Y f (CA ^Yl Q„»r. ; ,• l. r- U ->1 

5 

et on a /, (x, v) = .r - v ou /,(x,.y)=x/v 

*,(*)= 100.M ou ^ 2 (.r)=100.|log(|x|)| 

En particulier, la formuie BM-. = iqq. SA ^') est 
10 utilisee. SAjj) 

L'etape f met avantageusement en oeuvre la 
determination d'un indicateur 3M d'effet de bloc, selon 
la formuie suivante : 

SA d2 / SA-, 

15 BM = . x 100 

SA dl / SA X 

D'autres caracteristiques et avantages de 
1' invention apparaitront mieux a la lecture de la 
description qui va suivre, donnee a titre d'exemple non 
limitatif, en liaison avec les dessins ci-annexes, dans 
lesquels : 

la figure 1 est une representation 
schematique d'un dispositif pour la mise en oeuvre du 
procede selon la presence invention, 

la figure 2 represents une grille alignee 
pour ie calcul du parametre devaluation SA et une orille 
decalee pour le calcul du parametre SA., 

-a figure 3 represente 1' evolution du 
parametre BM au cours d 1 ur.e sequence filmee part iculiere, 

et les figures 4a a 4c representent 
respectivement devolution temporeile du parametre SA 
(4a) et la TCD appliquee a M echantiilons successifs de 
SA (4b et 4c) . 

Selon la figure 1, des images d' entree I. sont 
codees par un codeur COD et transmises par un systeme de 
transmission TRANS, par exempie un systeme de 
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transmission par satellite, puis decodees a la reception 
par un decodeur DEC pour donner une image de sortie 
ceccdee I s . Un circuit TCDE applique a titre d'exempie une 
transformee en cosinus discrete a 1 1 image d' entree I e et 
5 un circuit PARE determine les parametres Pe de 1 1 image 
d'entree I e qui seront utiles a devaluation a realises 
De meme, un circuit TCDS applique a titre d'exempie une 
transformee en cosinus discrete a 1 • image de sortie I s et 
un circuit PARS determine les "parametres P s de 1 1 image de 

10 sortie I s qui seront utiles a l'evaluation a realiser. 

Un circuit de comparaison COMP effectue des 
comparaisons entre les parametres d' evaluation ?e et Ps 
et produit en sortie des resultats devaluation R. Les 
parametres Pe sont transmis vers la sortie par ie systeme 

15 TRANS ou independamment de celui-ci. Alternativement , il 
est possible de faire "remonter" les parametres Ps vers 
1* entree afin qu'ils y soient compares avec les 
parametres Pe . 

On remarquera egalement qu'il est possible 

20 d'utiliser les coefficients TCD disponibles au niveau ' du 
codeur COD et/ou du decodeur DEC, mais alors l'evaluation 
ne tiendra pas compte des defauts du codeur COD et/ou du 
decodeur DEC. 

La transformee en cosinus discrete bi- 

25 directionnelle TCD fait correspondre a un bloc de N x N 
pixels (en general ,V « 8) de 1 ' espace pixel un bloc de N 
x iV coefficients du domaine transforme. Ces derniers 
expriment la representation f requentielle de l 1 amplitude 
de la luminance des pixels. 

30 La transformee TCD presente des coefficients 

F(u,v) dans iesquels : 

u,v = index de la frequence spatiale 
r.orizontale et verticale, 

x,y = index de position horizontale et 

35 verticale du pixel dans ie bloc, 
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- f(x,y) representee la luminance du pixel au 
point (x,y) dans 1" image, 

- F(u,v) represents le coefficient tranforme 

au DOinf (ll . vl rhi Hrvma ^ «- _ , i 

•• >— r - / —v- «wuiuin^ u j. a no i. u i. luc , 

5 " avec u,v,x,y variant de O a N-l (N = 8). 

Les parametres exploites pour representer les 
degradations contenues dans les images numeriques codees 
notamment en MPEG2 utilisent cette transformation. 

Une information sp'atiale est fournie a partir 
10 du domaine transforme par le calcul de l'activite 
spatiale globale SA de 1' image. Ce parametre informe 
precisement sur les details dans le signal image (l t ) 
indexe a 1 • instant t. (I t ) contient (HxW) blocs de pixels. 
La taille du bloc f m (x,y) est 8x8 pixels (voir fig.l). 

Dans un premier temps la transformer TCD par 
blocs est appliquee a 1 • image entierement, ou 
partiellement (sur une zone d'interet ou fenetre 
d'analyse). Le premier facteur representatif du contenu 
de !• image, l'activite spatiale, est obtenu a partir de 
1' image tranformee TCD (I c ). Un bloc image tres actif 
spatialement a une representation consequente dans le 
domaine transforme. Ses coefficients ont une amplitude et 
une distribution spectrale importantes . 

Deux parametres sont introduits, l'activite 
spatiale de bloc (bsa) et i'activice spatiale d ' image 
SA(I t ). Us represen-ent respect ivement l'activite 
spatiale intra-bloc et intra-image. 

Le contenu analyse est celui de chaque bloc de 
pixels f nm (x,y)de position a n,m dans 1 ' image (I,). 

L'indicateur (bsa n , ro ) pour le bloc n,m tient 
compte des composantes alternatives du bloc de 
coefficients transformers F n>m (u,v). 

On a : bsa„. m = ^(/_(.r, y)) = f £X{F n . m VJ))] "(F_(0,0)V 
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SA(I t ) est 1' expression du contenu global de 
1 1 image (ou de la zone d'anaiyse). Zl est calcule par une 
-oyenne ies activites spatiaies (bsa,^) des blocs 
concern.es . 



SA(I, ) = bsa^I, ) = — i— £ X M/-(*' v )) 



ou (H x W) est le nombre de blocs de pixels de 
10 l 1 image ou de la fenetre d'anaiyse. 

L'indicateur d'activite spatiale global SA(I t ) 
permet de qualifier la richesse de 1' image. Ce parametre 
est affecte par le codage a reduction de debit et par les 
erreurs de transmission. La comparaison du meme parametre 
15 calcule a 1 ■ entree et a la sortie d'un systeme video 
numerique contribue a 1 1 identification des degradations 
generees par ce systeme. Ce parametre est represenratif 
des degradations generees par le systeme, mais sans 
1 1 inconvenient d'une compensation entre les defauts. En 
effet, d'une part le calcul realise par blocs elimine les 
effets de bloc et d* autre part les pertes d' information 
et les faux contours se traduisent dans la tranformee par 
des effets de perte des hautes frequences, et ces deux 
phenomenes agissent done dans ie meme sens. 
25 L'evolution temporelle de SA(I t ) represence 

les evenements apparus dans la sequence video. Par 
ailleurs, ce parametre est affecie par le codage a 
reduction de debit et par les erreurs de transmission. La 
comparaison du meme parametre calcule a 1' entree et a la 
sortie d'un systeme video numerique est effectuee sur une 
image I t indexee a 1' instant t. Elle contribue a 
1 ' identification des degradations generees par ce 
systeme. La formulation mathematiaue utilisee oour la 



comparaison est . 
35 ! SA(l.) e 
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Le parametre definitif pour la mesure de la 
degradation liee a la perte d 1 informations est ie 
parameire iMLR : 

MLR M = Maximum.. [ log.rat io ( 1^ ) ] 
5 Cet indicateur est caicuie sur une sequence de 

M images. Cette option permet d 1 avoir une information 
reduite concernant une periode de temps. Elle permet 
d' avoir une note de qualite sur : T = M x 1/25 seconde, 
ou plus generaiement T = M/NI, NI designant le nombre 
10 d' images par seconde. 

La procedure de calcul de la perte 
d ' information se presence par exemple comme suit : 

a) selection d'une image a i'entree du 

systeme, 

b) calcul de SA(I t ) e pour cette image, 

c) indexation de cette image avec sa date et 
son parametre associe Sfc e , 

d) transmission de parametres de reference 
vers la sortie du systeme, 

f) calcul de SA(I C ) S pour cette image, 

g) calcul du log. ratio ( l z ) pour chaque image 



15 



20 



h) calcul de la noi:e ce qualite MLR pour 
chaque sequence de M images. 
25 On notera que le precede pourrait egalement 

etre mis en oeuvre en faisant "remonter" vers I'entree 
les valeurs de SA calculees pour chaque image de sortie. 

LJn avantage du procede est que chaque image a 
1' entree et a la sortie est caracterise par la valeur du 
30 parametre SA, ce qui signifie qu'il n'est besoin de 
transmettre en tant que parametre Pe qu'une seuie donnee 
par image pour caracteriser les degradations (autres que 
I'effet de bloc) par opposition a des procedes de 1 1 art 
anterieur qui necessitent de transmettre un grand nombre 
de donnees permettant d'effeccuer une comparaison par 



35 
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difference pixel a pixel entre les images reference et 
degradee . 

La cornparaison effectuee a partir du facteur 
MLR traduit une .tiesure de la qualite des images finales 

5 et par consequent informe sur les performances du systeme 
de codage mis en oeuvre lors de la transmission. 

Afin de valider la pertinence du facteur MLR, 
nous avons confronte les resultats obtenus aux notes 
obtenues par des panels de personnes representatives lors 

10 de tests subjectifs. Le fort coefficient de correlation 
obtenu entre 1' evaluation objective et subjective nous 
demontre I'efficacite d 1 un tel parametre. 

Dans les sequences video codees a faible 
debit, la distorsion par effet de bloc se manifeste par 

15 des discontinuity du contenu de 1 1 image aux limites des 
blocs adjacents. Cet artefact est visible sur les zones 
de faible ou de grande activite spatiale et aussi sur les 
zones de forte ou de faible luminosite. II est du a un 
codage inter-blocs independant et a une quantification 

20 grossiere du contenu de chaque bloc. Le coefficient 
continu F(0,0) de la transformee TCD represente la 
moyenne des valeurs de luminance d*un bloc. Le codage des 
coefficients continus de deux blocs, adjacents 
spatialement, a des valeurs sensiblement differentes est 

25 une des causes de cette degradation. Les autres 
coefficients transformes, qui representent les 
composantes alternatives, traduisent les differentes 
evolutions contenues dans le bloc. Le codage et la 
quantification grossiere de ces coefficients genere 

30 egalement 1' effet de bloc dans les zones des contours. 
Des erreurs dans le flux numerique influent sur la 
compensation de mouvement. Ceci se traduit par un effet 
au niveau des blocs ou des macro-blocs. 

Selon i' invention, un parametre de mesure de 

35 l 1 effet de bloc met en oeuvre une cornparaison de deux 
parametres calcules de chaque cote du systeme video. 
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L 1 indicareur propose permet c'evaluer 
1 ' importance de 1 1 information visueiie erronee apportee 
par la grille structuree de blocs 8x8 visible dans les 
images en sortie de chaxne. 
3 Les bordures des blocs 8x8 ou les 

discontinuities de contours sont considerees cornme une 
information perceptible et additionnelle dans 1' image. 
Sachant que le codage est realise de maniere independante 
pour les blocs adjacents, - les pertes d 1 information 

10 internes aux blocs generent des discont inuites de 
contenu, voire meme 1' apparition de nouveaux contours sur 
les frontieres. 

Un procede adapte a cette degradation selon 
1' invention repose sur le calcul de deux parametres : 

15 i'activite spatiale globale SA (voir supra) et I'activite 
spatiale globale decalee SA d (voir fig. 2). 

Le parametre SA reste non affecte par les 
effets de bloc, puisqu'il est realise a partir de ces 
memes blocs et de maniere alignee sur leur grille. Afin 

20 de souligner la contribution des contours et des defauts 
autour des blocs, on determine l'activite spatiale 
decalee SA d . Elle permet de mettre en exergue les 
artefacts aux frontieres. A partir de la meme definition 
de SA, SA d est calculee avec de preference un decalage de 

25 N/ 2 = 4 pixels hori zontalement et de N/ 2 = 4 pixels 
verticalement . On remarquera que l'activite spatiale 
decalee SA d peut etre calculee en effectuant un decalage 
de 1 a N-l pixels (et de preference N/ 2 pixels) seulement 
dans la direction horizontaie, ou seulement dans la 

30 direction verticale, de maniere a evaluer l'effet de bloc 
horizontal ou bien l'effet de bloc vertical. 

Les deux parametres sont issus d'une somme sur 
la zone d'interet ou sur 1 1 image entiere. Une comparaison 
entre les deux indicateurs met en valeur l'apport de 

35 1 1 information des contours au bord des blocs. Un nouveau 
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parametre, est done introduit, Sa ratlo dans 1' equation 
suivante. ^ (() 

5 S^ratio n'est pas seulemenr sensible aux contours ou 
discontinuites aux boras des blocs mais il l'est 
egaiement au contenu d'origine des blocs. II est par 
consequent necessaire, pour mesurer specif iquement la 
degradation introduite par- un systeme video, d'en 

10 extraire 1 1 information apportee par le contenu original 
de 1' image. 

L'indicateur definitif pour la mesure et la 
detection de l'effet de bloc introduit par un systeme 
video est calcule a partir des valeurs SA ratio obtenues a 
15 l f entree et a la sortie du-dit systeme. 

Ces valeurs sont considerees respectivement 
comme valeur du signal reference a 1' entree (SAratio e ) et 
valeur du signal degrade a la sortie (SAratio s ) . 
• L'indicateur de l'effet de bloc BM est defini par 
20 1' equation suivante. 

r>.„ 1,^ SAratio e 

BM SA = 1 00.— 

SAratio s 

Le procede de detection de i'effet de bloc 
25 peut etre mis en oeuvre par exemple comme suit : 

a) selection d'une image a 1 1 entree du 

systeme, 

b) calcul de SA(I t ) e , SA d {I.) e et SA rati0 (I t ) e 
pour cette image, 

30 c) indexation de cette image avec sa date et 

son parametre associe, 

d) transmission de donnees de reference vers 
la sortie du systeme, 

e) identification de la meme image indexee a 
35 la sortie du systeme, 
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f) calcul de SA(I t ) s# SA d (I,) s et SA ratlo (I t ) s 
pour la meme image a la sortie du systeme, 

g) calcul du BM SA pour chaque image I t , 

h) calcul de la moyenne de BM SA pour chaque 
5 sequence de M ' images (avec par exemple M' = M) . Cette 

note qualifie la sequence pour la degradation en 
question . 

La figure 3 demontre 1 ' importance decroissante 
de BM SA pour 4, 6 et respect ivement 9Mbits/s. Cette 

10 remarque est valable specialement dans la derniere partie 
de la sequence video. En effet, pendant cette periode le 
contenu spatial et temporel de la sequence choisie evolue 
brutalement. Ceci provoque une baisse de la cualite des 
images finales due a une reduction de debit. Le codage 

15 par quantification des coefficients TCD provoque 
1* apparition de 1' effet de bloc, et ce d'autant plus que 
le codage est grossier. 

Les degradations dans les sequences d' images 
numeriques codees ont des formes spatiales ou 

20 temporelles. En effet, des erreurs de codage ou de 
transmission peuvent affecter le contenu spatial : 
contours, objets et structures contenus dans le signal. 
Cependant, elles peuvent egalemenr entacher 1' evolution 
temporelle de la sequence video : saccades, gel d' images 

25 ou deplacements errones d ! objets. 

Afin d 1 identifier les degradations temporelles 
qui alterent les sequences video, nous avons defini un 
parametre supplementaire appele "Activite Temporelle" TA. 
II traduit les evolutions du contenu des sequences 

30 d* image dans le temps. II est sensible aux differents 
evenements qui peuvent survenir et affecter le signal. 

L'indicateur spatial SA(L) traduit le contenu 
de 1' image courante dans la sequence a un instant t. Les 
variations temporelles de SA(I t ) representent les 

35 evolutions temporelles du contenu de 1 1 image et done 
elles traduisent les changements survenus (voir fig. 4a). 
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Le signal mono-dimensionnel SA(I-) d'activite spatiale 
globale de 1 ' image est proport ionnel a la richesse du 
contenu. Dans le cas mono-dimensionnel, la TCD est une 
transformation qui fait correspondre a un ensemble de M 
points fix) une serie de M frequences. Eiie se definit 
par : 



F( U) = ^ -CM u)T/(x) cosi^— ; « . 



c(u) = \Nl sm = 0 
c(u) =1 si u * 0 

F(0)= S(0)S/(-v) 

V /V c-a 



Sur loute fenetre de M echantillons temporels 
15 du signal SA(I t ) est appliquee une transformee en cosinus 
discrete TCD. Cette nouvelle representation traduit ies 
frequences temporelles du signal et par consequent 
l'activite de la sequence d'image. La premiere raie F(O) 
est la composante continue, les autres raies F(u) avec 
20 u*o sont les composantes alternatives. 

La composante continue du spectre renseigne 
sur la valeur moyenne de la richesse du contenu sur M 
images . Les frequences alternatives informent sur ies 
variations temporelles du contenu. Seuie 1' amplitude des 
25 raies est prise en consideration. 

TCD T = TCD (SA(l t )) = ^F 0 (u) u = 0.. M-l 

Le spectre TC D T (voir fig. 4c) du signal SA 
30 (voir fig. 4b) informe sur les evolutions temporelles de 
la sequence video et par consequent traduit toute 
degradation modifiant ces evolutions, r (0) par definition 
represente la moyenne du signal SA(I.) pendant les M 
images de sortie considerees dans la sequence. Les 
35 amplitudes des composantes alternatives F(u) u ^ Q traduisent 
la maniere avec laquelle change le contenu de la video : 
variations rapides (clips video), modifications lentes 
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("travelling"), plans fixes, etc... Afin d' exploiter les 
proprietes de 1' analyse dans le domaine transformer un 
parametre TA a ete defini. II est: calcule pour un croupe 
de M images. ra " 

TA indique i ! amplitude des fluctuations de 
l'activite SA dans la sequence. II traduit tout evenement 
temporel de la sequence et parr consequent : les coupures 
10 et les gels d' images. TA est compare a un seuil donne S x 
etabli empiriquement , en dessous duquel on considere 
qu'on se trouve en presence d'une coupure ou d'un gel 
d ■ images . 

Le procede de detection des coupures et des 
15 gels d 1 image peut etre mis en oeuvre par exemple come 
suit : 

a) selection d'un groupe d' images, 

b) calcul de SA(I t ) pour chaque image, 

c) calcul de la TCD T et de TA, 

20 d) - si TA<s x et si F(0)<s 2 , alors le groupe 

d' images est une succession de la meme image uniforme (en 
pratique coupure noire) . 

- si TA<s x et si F(0)>S 2 , alors le groupe 
d 1 images est une succession de la derniere image non 

25 erronee gelee par le decodeur. 

- si TA>s i , retour a a. 

Le phenomene de gel d' images est une 
degradation caracterist ique qui apparait lors d'erreurs 
sur la chaine complete de distribution du signal de 
30 television. Le parametre d'activite temporelle TA est 
applicable dans ce cas de figure. En effet, il detecte 
cet evenement qui est lie a la strategie du decodeur. Le 
gel d' images prend deux formes : 

1. restitution des dernieres donnees images 
35 non corrompues, 
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2. visualisation d'une image unif.orme (souvent 
noire) , et ceci tant que les donnees regues sont 
inexploitables . 

La famille de codage MPEG introduit la notion 
5 de groupe d'images ("Group of Pictures GOP"). Cette 
structure de la sequence est liee a \ une syntaxe 
definissant les algorithmes pr^dictifs a ' reduction de 
debit qui sont appliques a chaque image (images codees 
sans prediction : Intra - I, images predictives 

10 unidirectionnelles P, images predictives Bidirectionelles 
B) . La particularity de chaque image influe sur I 1 aspect 
de 1' impact des differentes erreurs qui I'affectent. 

L'indicateur defini ici traduit les evenements 
temporels qui surviennent durant la sequence video. II 

15 permet par consequent d 1 identifier et de comptabiliser 
les ruptures temporelles des sequences. Par la meme 
occasion il represente un parametre de qualification de 
la qualite du service. 

L 1 invention ne se limite pas aux modes de 

20 realisation decrits, En particulier , l'effet de bloc 
pourrait etre evalue de naniere intrinseque sur 1 ' image 
de sortie. Le procede est alors caracterise en ce qu'il 
comporte les etapes suivantes : 

di) determiner I'activite spatiaie globale SA d2 

25 de l 1 image de sortie dans une fenetre d* analyse decalee 
presentant des blocs de points image ou pixels qui sont 
decales d ! au moins un pixel dans la direction des lignes 
de 1* image et/ou dans une direction perpendiculaire a la 
direction des lignes de 1 1 image par rapport auxdits blocs 

30 de pixels de ladi.e fenetre d' analyse, cette 
determination mettant en oeuvre les sous-etapes suivantes: 
i) appliquer ladite transformee par blocs sur 
chaque bloc decale (n,m) dudit ensemble de blocs de 
pixels de la fenetre d 1 analyse decalee pour determiner 

35 les coefficients transformers F d n,m (i/j), 
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ii) determiner a partir des coefficients 
transformes F d n,m (i,j) des blocs decales, l'activite 
spatiale bsa d de chaque bloc dudit ensemble de blocs 
decales, 

iii) determiner l'activite spatiale globale 
SA^ de 1' ensemble des blocs decales constituent la 
fenetre d' analyse decalee, 

ei) comparer l'activite spatiale globale SA d2 
de 1' image de sortie dans la fenetre d • analyse decalee et 
l'activite spatiale globale SA 2 de 1 ' image de sortie dans 
la fenetre d' analyse. 

L' invention peut etre mise en oeuvre a l'aide 
de tranformees par blocs lineaires et orthogonales telles 
que la tranformee en cosinus discrete TCD, la transformee 
15 en sinus discrete, la transformee de Fourier, la 
tranformee de Haar, la transformee d'Hadamar, la 
transformee de Slant ou bien encore une transformee en 
ondelettes (DWT) dont les coefficients ont ete 
reorganises pour obtenir les blocs transformes. 
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REVENDICATIONS 

1. Proceae d' evaluation de la degradation 
d'une image video codee par blocs de points images ou 
pixels, cette degradation etant generee par un systeme de 
5 codage et/ou de stockage et/ou de transmission, 
produisant une image de sortie a partir d'une image 
d'entree, caracterise en ce qu'il comporte les etapes 
suivantes : 

a) selectionner une image d' entree et 
10 determiner une activite spatiale (SA) de 1 1 image d' entree 

dans une fenetre d' analyse representant au moins une 
partie de 1 ' image qui presente un ensemble de dits blocs 
de points images ou pixels, cette determination mettant 
en oeuvre les sous-etapes suivantes : 
15 i) determiner pour chaque bloc (n,m) du dit 

ensemble de blocs de pixels, les coefficients transformes 
Fn,m (i,y) par une transformee par blocs selon les 
dits blocs 

ii ) determiner a partir des coefficients 
20 transformes Fn,m (i,y) l'activite spatiale bsa de chaque 

bloc dudit ensemble de blocs de pixels 

iii ) determiner a partir de 1' activite 
spatiale de bsa de chaque bloc 1' activite spatiale 
globaie (SA X ) de i 1 ensemble des blocs de pixels 

25 constituant la fenetre d' analyse 

b) selectionner 1 1 image de sortie 
correspondant a 1 1 image d 1 entree et determiner ladite 
activite spatiale globaie (SA 2 ) de 1 ' image de sortie dans 
ladite fenetre d 1 analyse, en mettant en oeuvre pour 

30 1 1 image de sortie les sous-etapes a i) a a iii) 

c) comparer 1' activite spatiale globaie (SA 2 ) 
de 1 1 image de sortie dans la fenetre d 1 analyse et 
1' activite spatiale globaie (SAO de 1 * image d' entree dans 
la fenetre d'analyse. 
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2. Procede selon la revendicaiion 1 
caracterise en ce que l'activite spatiale bsa c'un bloc 
(n,m) est ceterminee par la formule suivante : 

5 bsa njn Jii{kU,j).F njn {i,j)) P ) -(^(0.0).F njJ ,(0.0)) P 

V<=0/=0 J 

ou k(i,j) est un coefficient constant de 'la ponderation 
pour la composante i,j utilisee et p une constante. 

3. Procede selon la revendication 2 
10 caracterise en ce que l'activite spatiale SA de 

1' ensemble des blocs peut etre alors determinee par la 
formule suivante : 

15 

(HxW) representant le nombre de blocs de 
pixels dans la fenetre d f analyse. 

4. Procede selon la revendication 3 
caracterise en ce que ladite comparaison (etape c) est 

20 effectuee a I'aide d'un parametre LR defini de la maniere 
suivante : 

LR = gj [fi (SAi/SA 2 ) ] avec i,je{l,2} 

avec fi(x,y)=x-y ou f2(x,y)=x/y 

et g:(x)-100lx| ou g2(x)=100 I log ( I x I )| 

25 

5. Procede selon la revendication 4 ■ 
caracterise en ce qu'il met en ceuvre le calcul du 
parametre LR sur une sequence de M images d' entree et de 
M images de sorties correspondants, et en ce qu'il met en 

30 oeuvre le calcul d'un parametre devaluation ou note de 
qualite MLR defini de la maniere suivante : 
MLR = Maximum* (LR) 

6. Procede selon une des revendications 
precedentes caracterise en ce qu f il comporte la 

35 determination d 1 un indicateur d'activite temporelle TA 
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pour un groupe de M images, determine de ia maniere 
suivante : r-jj-j 

^sa( u )^ pour u variant de 0 a M-l, designant 
les M coefficients d'une transformee par blocs appliquee 
a M echantillons temporels de l'activite spatiale SA. 

7. Procede selon la revendication 6 
caracterise en ce qu'il compofte une etape de comparaison 
de l'indicateur d'activite temporelle TA a un seuil donne 
S x , de maniere a determiner la presence d'une coupure ou 
d 1 un gel d' image. 

8. Procede selon une des revendications 
precedentes caracterise en ce qu'il comporte, en vue de 
la determination de maniere intrinseque de l'effet de 
bloc, les etapes suivantes : 

di) determiner l'activite spatiale globale SA d2 
de 1 * image de sortie dans une fenetre d' analyse decalee 
presentant des blocs de points image ou pixels qui sont 
decales d'au moins un pixel dans la direction des lighes 
de 1 1 image et/ou dans une direction perpendiculaire a la 
direction des lignes de 1' image par rapport auxdits blocs 
de pixels de ladite fenetre d 1 analyse, cette 
determination mettant en oeuvre les sous-etapes suivantes 

i) appliquer ladite transformee par blocs sur 
chaque bloc decale (n,m) audit ensemble de blocs de 
pixels de la fenetre d* analyse decalee pour determiner 
les coeeficients transformes F d n,m (i,j), 

ii) determiner a partir des coefficients 
transformes F d n,m (i,j) des blocs decales, l'activite 
spatiale bsa d de chaque bloc dudit ensemble de blocs 
decales , 

iii) determiner l'activite spatiale globale 
SA d2 de l 1 ensemble des blocs decales constituant la 
fenetre d ' analyse decalee, 
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ei) comparer I'activite spatiale globale SA d2 
de 1 1 image ce sortie dans la fenetre d'analyse decaiee et 
I'activite spatiale globale SA 2 de 1 1 image de sortie dans 
la fenetre d'analyse. 

5 9. Procede selon une des revendications 

precedentes caracterise en ce qu ! il comporte, en vue de 
la determination de I'effet de bloc, les etapes, suivantes 

dz) determiner 1' activity spatiale globale SA dl 
de 1 ' image d'entree dans une" fenetre d' analyse decaiee 

10 presentant des blocs de points image a pixels qui sont 
decales d'au moins un pixel dans la direction des lignes 
de 1' image er/ou dans une direction perpendiculaire a la 
direction des lignes de 1 ' image par rapport auxdits blocs 
de pixels de ladite fenetre d'analyse, cette 

15 determination mettant en oeuvre les sous-etapes suivantes 

i) appliquer une transformee par blocs sur 
chaque bloc decale (n,m) dudit ensemble de blocs de 
pixels de la fenetre d'analyse decaiee pour determiner 
les coeeficients transformes F d n,m (i,y), 

20 ii) determiner a partir des coefficients 

transformes F d n,m (i,y) des blocs decales, 1'act.ivite 
spatiale bsa d de chaque bloc dudit ensemble de blocs 
decales, 

iii) determiner I'activire spatiale globale 
25 SA dl de 1' ensemble des blocs decales constituant la 
fenetre d'analyse decaiee, 

e) determiner I'activite spatiale SA d2 de 
1' image de sortie dans ladite fenetre d'analyse decaiee 
en mettant en oeuvre les sous-etapes d 2 i) ad iii) 
30 f) comparer d'une part I'activite spatiale 

SA d2 de 1' image de sortie dans la fenetre d'analyse 
decaiee et I'activite spatiale SA 2 de 1 ' image de sortie 
dans la fenetre d'analyse, pour evaluer I'effet de bloc 
dans 1 ' image de sortie et d' autre part I'activite 
35 spatiale SA dl de 1 ' image d'entree dans la fenetre 
d'analyse decaiee et I'activite spatiale SA l de 1 1 image 
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d'entree dans la fenetre d'analyse pour evaluer l'effet 
de bloc dans 1' image d 1 entree. 

10. Procede selon la revendicat ion 9 
caracterise en ce que I'etape f met en oeuvre la 
determination d'un indicateur BM d'effet de bloc, selon 
la f ormule suivante : 

BM = x 100 

SA dl / SA X 

11. Procede selon une des revendications 
precedentes caracterise en ce que ladite transformee par 
blocs est une transformee lineaire et orthogonale, 
notamment une transformee en cosinus discrete TCD ou une 
transformee en sinus discrete ou une transformee de 
Fourier ou de Haar ou de Hadamar ou de Slant. 

12. Procede selon une des revendications 1 a 
10 caracterise en ce que ladite transformee par blocs est 
une transformee en ondelettes (DWT) dont les coefficients 
ont ete reorganises pour obtenir des blocs transformes. 
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